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Abstract 

Squark mixing plays a large role in the phenomenology of the minimal supersymmetric standard 
model, determining the mass of the lightest Higgs boson and the electroweak interactions of the 
squarks themselves. We examine how mixing may be investigated in high energy e + e~ reactions, 
both at LEP-II and the proposed linear collider. In particular, off-diagonal production of one 
lighter and one heavier squark allows one to measure the squark mixing angle, and would allow 
one to test the mass relations for the light Higgs boson. In some cases off-diagonal production may 
provide the best prospects to discover supersymmetry. In the context of the light bottom squark 
scenario, we show that existing data from LEP-II should show definitive evidence for the heavier 
bottom squark provided that its mass < 120 GeV. 
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I. INTRODUCTION 



Scalar quarks, the supersymmetric partners of ordinary colored fermions, are an impor- 
tant ingredient in any theory which combines the standard model (SM) with supersymmetry 
(SUSY). These scalar partners of the quarks play a key role in the softening of quadratic 
divergences in the Higgs boson mass parameter induced by loops of quarks. This reduced 
sensitivity to the ultra-violet physics is the single most attractive feature of the minimal 
supersymmetric standard model (MSSM). Phenomenologically, scalar quarks (squarks) pro- 
vide interesting signatures at colliders. They are color triplets, and so may be produced 
copiously in hadronic reactions. They also couple electroweakly, and thus may have spec- 
tacular decays involving leptons and weak bosons. In fact, as quarks in the SM are the fields 
which 'bridge' the strong and electroweak interactions, so squarks are the bridge between the 
supersymmetric analogues of the gauge sectors. As spin-0 fields, the angular distributions 
of their production and decay are often quite distinct from their spin-1/2 counterparts. 

Each massive quark of the SM results from the marriage (through electroweak symmetry- 
breaking) of two Weyl fermions with very different electroweak characteristics. Each Dirac 
fermion is accompanied by two complex scalars, also with distinct electroweak interactions. 
The supersymmetric versions of the Yukawa interactions between the Higgs boson and the 
fermions mix this pair of scalars, leading to two mass eigenstates denoted q\ and q\. The 
amount of mixing, represented in terms of a mixing angle 8 q , is typically proportional to the 
associated fermion's mass, and the mixing is thus presumably largest for the third generation 
squarks. 

Squark mixing plays an interesting role in the phenomenology of the MSSM. Perhaps 
foremost is the fact that the lightest MSSM Higgs boson, the one primarily responsible for 
electroweak symmetry-breaking, has a mass at tree level which is less than the Z boson 
mass, considerably smaller than 114.5 GeV, the experimental bound from the CERN Large 
Electron Positron (LEP) facility [1]. This expectation would exclude the MSSM if it were 
not for the fact that the Higgs boson mass receives large radiative corrections which can lift 
the mass beyond the reach of the exclusion limits, to masses as large as about 135 GeV. 
These necessary corrections are provided largely by the top squarks and are enhanced by 
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the squark mixing [2], 

*(!JHM)]- 

In this expression, m t is the top quark mass, Ms is a common scale of SUSY-breaking 
which describes the overall magnitude of the squark masses, and x t = (A t — fi cot P) /Ms 
characterizes the off-diagonal entries in the squark mass matrix (and thus the amount of 
mixing; x t oc sin6> t cos 6^), with the largest Higgs boson masses realized for x t ~ VQ- Once 
the lightest Higgs boson and top squarks are discovered, a key test that the Higgs boson 
properties follow the MSSM paradigm will require careful measurement of the top squark 
masses and their mixing angle. When this dominant correction is understood, one begins 
to indirectly probe the remainder of the MSSM through the sub-dominant one-loop and 
dominant two- loop corrections to the Higgs boson mass [3] . 

The mixing angle determines the couplings of the squark mass eigenstates to the W and 
Z bosons. It can be determined efficiently if these couplings are measured. For instance, 
the rate of the reaction e + e~ — > q*q2, production of one lighter and one heavier squark, is 
proportional to the amount of mixing, and a measurement of the rate is an excellent way 
to establish the value of the mixing angle. Off-diagonal production provides an interesting 
complement to other proposals [4] to measure the mixing angle at an electron-positron 
linear collider [5]. These measurements typically involve the cross section for top squark 
tit\ pair production with various combinations of polarized e + and e~ beams. Off-diagonal 
production provides an important cross-check of these other methods, and it does not rely 
on any particular polarization of the incoming beams. Even if the heavier squark mass is 
larger than half of the collider energy, precluding pair production of the heavier squarks, 
this method can still succeed and also allow one to measure the mass of the heavier squark. 

This line of reasoning is indicative of a general feature of squark phenomenology. While 
hadron colliders generally offer the best prospects to discover squarks, due to copious produc- 
tion through the strong interactions, it is somewhat problematic to determine parameters 
in hadronic reactions. In hadron collisions, the primary means to study the electroweak 
properties of squarks is to examine squark decays. This task is challenging for a number of 
reasons. In many popular models of SUSY-breaking [6], the squark decay is dominated by a 
single channel (often into a quark and the lightest neutralino) dictated purely by kinematics, 
and not by couplings. The branching ratios of various decay modes are relatively insensi- 
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tive to the squark mixing angle. Furthermore, the squark decay width, the quantity that 
actually depends on the coupling strength, is in almost all cases well below the experimental 
resolution and thus unmeasurable at hadron colliders. 

This discussion illustrates the importance of squark mixing on the determination of 
bounds on squark masses from data at lepton colliders. A particular example is furnished by 
the scenario of light bottom squarks and light gluinos proposed in Ref. [7] . In this scenario, 
the excess rate of bottom quark production at hadron colliders is explained by postulating a 
tree-level contribution from production of light gluinos which decay into a bottom quark and 
a bottom squark. Data indicate that the lifetime of the hypothesized light bottom squark 
must be less than 1 nanosecond [8]; in typical collider detectors, it does not have a signifi- 
cant missing energy signature nor does it produce tracks characteristic of heavy long-lived 
objects. The light bottom squark b\ is assumed to decay hadronically, via R-parity violation, 
without a visible flavor tag necessarily. Its signals are extremely difficult to extract from 
backgrounds [7, 9]. Furthermore, a light bottom squark evades LEP-I data if one uses the 
freedom to select a mixing angle which renders its coupling to the Z boson tiny [10]. This 
essential requirement implies a non-zero mixing angle (sin 2 9b ~ 1/6 ), and, therefore, the 
off-diagonal Z-b\-b 2 coupling must be non-zero. Here b 2 denotes the heavier of two bottom 
squarks. One of the most promising (and potentially clear) signals of this scenario is fur- 
nished by e + e~ — > Z* — > b x b 2 . References to other work on the phenomenology of light 
bottom squarks and light gluinos may be found in Refs. [9, 11]. 

In this article we consider off-diagonal squark production at e + e~ colliders. We compute 
tree-level cross sections and discuss likely decay modes. In Sec. II, we review the squark mass 
matrix and electroweak interactions. In Sec. Ill we compute the production cross sections 
for both bottom squarks and top squarks, showing the dependence on the masses and mixing 
angles. In Sec. IV we apply our results to the light gluino and bottom squark scenario at 
LEP-II, estimating for the first time the discovery potential of the heavier bottom squark 
to be greater than 5 standard deviations (5er) provided m~ b2 < 120 GeV. Alternately, the 
LEP-II data can exclude masses smaller than 130 GeV, if no signal is observed. We reserve 
Sec. V for conclusions. 
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II. SQUARK MASSES AND ELECTROWEAK INTERACTIONS 



In this section we briefly review the squark mass matrices, mixing angles, and electroweak 
interactions. The squark mass matrices are 

\m t {A t - ^ cot (3) M 2 R +m 2 t +D R ' 1 

for top squarks, and 
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Ml + ml + DL ^-" tan « 1,^, (3) 
m b (A b - fitanp) M 2 R +m 2 + D R 

for bottom squarks. Parameters M L and M qR are the SUSY-breaking masses for left- and 
right-handed squarks, A t and A b are the (SUSY-breaking) trilinear interactions with the 
Higgs field, n is the Higgsino mass parameter, tan/5 is the ratio of Higgs boson vacuum 
expectation values (VEVs), m b and m t are the bottom and top quark masses, and and 
and D R are D-terms, D L = m| cos2/3(T 3 — Qf sin 2 9w), and D R = Qffn 2 z sin 2 6 W cos 2(3, with 
Qf the quark charge, and T 3 its weak isospin. For simplicity, we neglect the possibilities of 
phases in the soft-breaking parameters; their inclusion is straight-forward. 

Diagonalization of the matrices determines the squark mass eigenstates characterized by 
mixing angles 9 q . The physical squarks are a mixture of the scalar partners of the left- and 
right-chiral quarks. The mass eigenstates are two complex scalars (q± and ^2), expressed in 
terms of left-handed (L) and right-handed (R) squarks, and q R , as 

|<?i) = sm9 q \q L ) +cos9 q \q R ), 

\q 2 ) = C os9 q \q L ) - sin9 q \q R ), (4) 

where our convention is that q\ is the lighter of the two mass eigenstates. These angles may 
be expressed as 

. 2m t (A t -//cot/3) 

sm2^ = k ^ , (5a) 

m~ — mj 

. 2m b {A b -//tan/3) 

sm26> 6 = 5 5 . (5b) 

mf — m{ 
01 02 

It is worth emphasizing that the mixing terms in Eq. (3) are proportional to the quark 
masses, and thus the mixing is presumably largest for the third generation squarks. For this 
reason, we focus on scalar top and bottom quarks in the discussion below. 



The electroweak interactions of the squarks are determined by the relative admixture 
of left- and right-chiral squark in the mass eigenstate. The coupling to the SU(2) L gauge 
bosons is only through the left-chiral component, whereas coupling to the U(1) Y boson 
is non-zero for both. After electroweak symmetry breaking (EWSB), the photon remains 
massless, and its gauge invariance is linearly realized. All squarks (of a given electric charge) 
couple equally to the photon with coupling strength given by Qfe. The Z boson couplings, 
on the other hand, are sensitive to the mixing angles, 

sin aw cos u\y 

g 

# 12 = ~ ~a p3sin0 9 cos0,] , 

sm Ow cos Ow 

g 22 = 6 [T 3 cos 2 9 q - Q f sin 2 9 W ] , (6) 

sm a\y cos u\y 

where gij refers to the coupling of the Z boson with ^ and cjj. As mentioned above, the 
coupling to the lighter squarks may be tuned to vanish for sin 2 9 b ~ 1/6 (for bottom squarks) 
and sin 2 4 ~ 1/3 (for top squarks) [10], but in these limits, the off-diagonal couplings, and 
the heavy- heavy couplings are non-zero. 
The couplings to the W boson are 
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where now refers to the coupling of W to U and bj. Squark couplings to quarks and 
either charginos or neutralinos x° are straightforward, but somewhat more complicated 
by the mixing angles associated with the x* 1 and x° mass eigenstates. 



III. PRODUCTION CROSS SECTIONS 



In this section, we examine the production of the off-diagonal pair of squarks q\ and q\ in 
the electron-positron annihilation process e + e~ — > qiq%, illustrated in Fig. 1. The conjugate 
process e + e~ — > q{(j2 has the same cross section. Each reaction has a single Feynman diagram 
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FIG. 1: Feynman diagram for the process e + e — > qiq 2 . 

in which a Z-boson in exchanged in the s-channel. The unbroken gauge invariance of QED 
forbids the photon from contributing to off-diagonal squark production. For our purposes, 
it is enough to consider the tree-level production rates. Initial state radiation, and Yukawa 
and SUSY-QCD one-loop corrections are computed in Ref. [12] and can be typically as large 
as ±15% for some regions of parameter space. 



A. e+e — > q 1 q% 

The amplitude for the process 

e-(fci) + e + (k 2 ) -> Z* ^b 1 (p 1 ,m 1 ) + b*(p2,m 2 ) (8) 

is expressed as 

M = 9uv(k 2 ) lfl (g R P R + g L P L )u(h) (9) 

z 

where Pl — (1 — 7s)/2 and Pr = (1 + 7s)/2. Specified are the four-momenta k\ and k 2 of 
the incident e~ and e + , and pi and p 2 of the final squarks. The lepton couplings to the Z 
are 

9l = . a 6 2 (-1 + sin2 9 w) 5 (10a) 
sin 9 W cos 9 W \ 2 / 

9R = ■ n 6 n (sm 2 ^) • (10b) 
sin V w cos U w 

Taking the absolute square of the amplitude, summing over final spins and colors, and 

averaging over the initial spins, we obtain the differential cross section 

dv _ r l2= Sgj 2 (g 2 L + g R ) r 3 sm 2 6* 

dcosO* 32irs^ 1 1 128tts (1-Mf/s) 2 ' 

where 



Angle 9* is the scattering angle of q~i in the e + e~ center-of-mass frame, and p* is its three- 
momentum. As expected for production of scalar particles, the energy dependence of the 
cross section is influenced by a P-wave threshold factor oc |p*| 3 , and the angular distribution 
varies as sin 2 9* . 

The rate integrated over the scattering angle is 

a = iM±M — ^_ (13) 

32tts (l-M 2 /s) 2 ^ ' 

in which the dependence on the squark mixing angle 9 q is manifest in the proportionality 
to sin 2 29 q (in the factor of g\ 2 ) in Eq. (11). For reference, note that production of a q± or a 
q~2 pair proceeds through both photon and Z exchange, and is not proportional to sin 2 29 q , 
containing terms both sensitive and insensitive to the mixing angle. 

B. Rates at a Linear Collider 

As discussed above, a key measurement at a future linear collider would be to verify 
the Higgs boson mass dependence on the supersymmetry-breaking parameters. This test 
would demonstrate that the MSSM is the effective theory at the weak scale, as opposed to 
some more general supersymmetric extension. As in most of the MSSM parameter space, 
the dominant corrections to the Higgs boson mass are from the top squarks, it is their 
electroweak properties that are most relevant. We envision (for illustrative purposes) a 
situation in which the Higgs boson has been discovered at the Large Hadron Collider through 
some combination of production and decay channels (see Ref. [13]), and the squarks have 
been produced through the strong interaction, dominantly gg — > t*ti and gg — > t^2- The 
squark masses and dominant decay channels are likely to be known, but the mixing angle 
(which plays no role in the tree-level production through the strong force) must be measured 
at a linear collider. 

We consider, for reference, the light top squark to have a mass of 150 GeV, somewhat 
above the Fermilab Run I bounds [14] , although the bounds themselves are sensitive to the 
details of how the top squark decays. In Fig. 2, we show the cross sections for e + e~ — > t*t2 
as a function of the heavier top squark mass, for a reference mixing angle of sin 2 9 t = 1/2. 
We choose three center-of-mass energies: y/s = 500, 800, and 3000 GeV. The rate is doubled 
if the charge conjugate process is included. We see that rates are typically of order a few 
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FIG. 2: Cross section for the process e + e~ — ► t\i\ at center-of-mass energies 500, 800, and 3000 
GeV, as a function of the mass of the heavier top squark and for mixing angle sin 2 Ot = 1/2. The 
mass of the lighter top squark has been fixed to 150 GeV. 



femtobarns (fb) for top squark masses within the range allowed by kinematics. A linear 
collider with hundreds of inverse fb of data could expect to produce (before cuts) hundreds 
of events, and the cross section could be measured at the few per cent level, provided 
experimental efficiencies are not extremely small and backgrounds not prohibitively large. 
Such questions must be answered in the context of specific top squark decay signatures and 
are not addressed in this work. 



IV. LIGHT BOTTOM SQUARKS AND LEP-II 

As our second example, we consider bottom squark production and adopt parameters 
suggested in the light bottom squark scenario of Ref. [7]. This scenario postulates that the 
excess rate of bottom quark production at hadron colliders arises from pair production of 



9 



gluinos with masses on the order of 15 GeV. The gluinos subsequently decay to bottom 
quarks and the light bottom squarks, with masses of order the bottom mass. In order for 
such light scalar bottom quarks to be consistent with Z-pole data, the light bottom squarks 
must decouple from the Z, implying a non-trivial mixing angle sin 2 9b ~ 1/6. Thus, the 
off-diagonal coupling to the Z boson is necessarily non-zero. 

The viability of the light bottom squark scenario has been questioned on the grounds 
that the heavier b 2 should have been detected at LEP-II. The argument is based on the 
evaluation of SUSY-QCD corrections to the Zbb vertex in Ref. [15] within the context of the 
light bottom squark and light gluino scenario. These loop corrections contribute negatively 
to Rb, the ratio of the width for Z — > bb to the total hadronic width, and they increase in 
magnitude with the mass of b 2 . To maintain consistency with data, the authors of Ref. [15] 
argue that the mass of b 2 must be less than 125 (195) GeV at the 2a (3a) level. In an 
extension of this analysis, Cho claims that b 2 must be lighter than 180 GeV at the 5a 
level [16]. A b 2 in the mass range < 200 GeV could have been produced in association 
with a light b\ at LEP-II energies, and since no claim of observation has been made, the 
authors of these studies suggest that LEP data disfavor the light bottom squark scenario. A 
heavier b 2 200 GeV) is allowed if CP- violating phases are present [17]. Real decays such 
as Z — > b\bg + b\bg contribute positively to Rb [18] and soften these bounds to 160 (290) 
GeV at the 2a (3a) level [19]. We remark that experimental searches for SUSY particles are 
model-dependent and a search of LEP-II data for a b 2 in the light bottom squark scenario 
has not yet been undertaken. The cross sections and discussion of decay modes in this paper 
may help to motivate such a search. 

We begin with the predicted cross sections and event rates for production of bib 2 pairs 
at energies explored at the CERN LEP collider. For the large heavy bottom squark masses 
that we consider, LEP-II is unable to pair-produce heavy bottom squarks, and off-diagonal 
production is the only viable option. We then discuss decay of b 2 , the heavier of the two 
bottom squarks. We consider the dominant decay mode b 2 — > bg, and we evaluate the total 
width for this decay as a function of m 2 , the mass of b 2 . Subsequently, taking gluino decay 
into account, we present and evaluate the amplitude for the full three-body decays b 2 — > bbb\ 
and b* 2 — > bbb\. The Majorana nature of the gluino permits final states in which there can be 
bottom quarks of the same sign (i.e., bb or bb) as well as the bb configurations expected in 
SM situations. The overall process, e + e~ — > b\b* 2 , followed by b 2 decay leads to a four-parton 
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120 150 180 m 2 (GeV) 

FIG. 3: Cross section for the process e + e~ — ► bib^ at center-of-mass energies 130, 161, 189, and 
207 GeV, as a function of the mass rri2 of the heavier bottom squark. The mass of the lighter 
bottom squark is m\ = 3.5 GeV, and sin 2 9b = 1/6. 

final state. 

A. Cross Sections and Event Rates 

Selecting center-of-mass energies spanning those at which data were accumulated at the 
CERN LEP-II facility, we show the cross section for b]b* 2 production as a function of the 
mass 777,2 in Fig. 3. In this illustrative calculation, the mass mi of the lighter bottom squark 
is mi = 3.5 GeV, and sin 2 9b = (2/3) sin 2 9w — 1/6. Focusing on the energy dependence 
at m 2 = 100 GeV, we notice that the cross section grows with center-of-mass energy y^i 
from 130 to 161 GeV and then falls as energy increases. This behavior may be traced to 
the combined influences of the |p*| 3 threshold suppression and the usual 1/s dependence at 
large s. 

Multiplying by the accumulated integrated luminosities per experiment [20] at LEP-II, 
we use our cross sections to compute the predicted number of events produced as a function 
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120 150 180 m 2 (GeV) 

FIG. 4: Number of events from e + e~ — > 6162 at center-of-mass energies 189, 196, and 207 GeV as 
a function of the mass of the heavier bottom squark. The mass of the lighter bottom squark is 
m\ = 3.5 GeV, and sin 2 8^ = 1/6. 

of rri2- These results are shown in Fig. 4. Below LEP-II center-of-mass energy 189 GeV, the 
integrated luminosities were too small to have produced an appreciable sample of events for 
the process of interest to us. In order to translate the event rates in Fig. 4 into limits on 
the observability of 62 ; we must discuss likely decay modes, experimental efficiencies, and 
backgrounds. Decays are discussed in the next subsection. Here we remark simply that if 
at least 5 events are deemed necessary, the raw event rates in Fig. 4 suggest that a bottom 
squark with mass greater than 140 GeV will have escaped detection at LEP-II. We present 
more rigorous estimates below. 

B. 62 Decay 

In a scenario in which the gluino g is lighter than 6 2; the most likely decay process is 
62 — > bg. As derived below, both the 62 and g widths are narrow compared to their masses, 
and thus the description of the three-body decay 6 2 — > bbb\ in two sequential steps, 62 — ► gb 
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followed by g — > bb\ is an accurate one. The width for b 2 — > gb, in the mj = limit, is 

In this expression, denotes the gluino mass. In the limit m 2 ^> m h) m gi the width 
grows linearly with m 2 , as expected. To estimate the magnitude of r(6 2 — *• <?&), we adopt a 
gluino mass within the range obtained in the light gluino and light bottom squark scenario: 
12 < m g < 16 GeV. The full width is more than an order of magnitude smaller than the mass 
m 2 , as expected since the relative size is controlled by a s (m 2 ) ~ 0.1. For example, choosing 
m 2 = 150 GeV and trig = 15 GeV, we find 1^ ~ 10 GeV. In our subsequent treatment of the 
process e + e" — > 6163, with b 2 — > bbb \ or b 2 — > bbbi, we are justified in adopting the narrow 
width approximation for b 2 , factorizing the production and decay. 

In the light gluino and light bottom squark scenario, the gluino decays with 100% branch- 
ing fraction into a bottom quark and and a light bottom squark. However, since the gluino 
is Majorana in nature, it may decay into either a bottom quark or a bottom anti-quark: 
g — > bb\ or g — > bb\. As an intermediate step in our full calculation of the width for the 
three-body decay of 6 2 , we first evaluate the width for on-shell gluino decay. The decay 
width for the two-body subprocess g — > bb is 

T g = T{g -> 66;) + T(g - 66^ = m„ (15) 

where the small bottom and light bottom squark masses are neglected. For m g = 15 GeV, 
we find T g = 0.6 GeV. We note that corrections to this expression from the finite bottom 
quark and bottom squark masses are generally not negligible, and introduce a dependence 
on the squark mixing angle. However, in all cases T g m g , and these corrections have little 
effect on the heavy bottom squark width or the distributions of the (highly boosted) gluino 
decay products from 6 2 decays. 

C. b* bb b* and b* -► bb 6i 

In this subsection, we address the full three-body decay subprocesses 6 2 — > bbb\ and 
b* 2 — > 666i. The relevant Feynman diagram in the case of opposite sign 'OS' production (66) 
is shown in Fig. 5, and the two diagrams for like-sign 'LS' production in Fig. 6. For OS and 
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TABLE I: The decay widths IYs and Tos hi GeV obtained with m g =15 GeV, sin 2 Ob = 1/6, and 
a s {m 2 ) =0.116, 0.112, 0.109 and 0.108 for m 2 =100, 125, and 150, and 175 GeV, respectively. Also 
shown is the total I\ . 

m 2 = 100 GeV 125 GeV 150 GeV 175 GeV 

T LS 3.8 GeV 4.6 GeV 5.4 GeV 6.2 GeV 

r OS 3.6 GeV 4.5 GeV 5.4 GeV 6.2 GeV 

7.4 GeV 9.1 GeV 10.8 GeV 12.4 GeV 



b 




g 



b 



FIG. 5: Feynman diagram for the decay b\ — > bbb\. 
LS decay, our kinematic notation is 

OS : b*(s 2 )^b(p 1 )+b(p 2 )+bl(p i ); (16a) 
LS : ~b*(s 2 ) -.6(pi) +b(p 2 ) +6i(pg), (16b) 

where the quantities in parenthesis are 4-momenta. In evaluating the amplitudes for these 
decays, we must contend with the fact that the gluino goes onto its mass-shell within the 
physical region. To handle this singularity, we resum a class of contributions to the imaginary 
part of the gluino 2-point function to all orders, replacing the gluino propagator by a Breit- 
Wigner form, p 2 — m~ — > p 2 — m~ + irrigYg, and we use the expressions for Y g above. Explicit 
expressions for the matrix elements and decay widths are presented in the Appendix. 

b 

b 

b; 





(a) (b) 
FIG. 6: Feynman diagrams for the decay b^ — > bbb\. 
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Defining 



Tls = r(65->6i+66), (17a) 

r os = rm^h + bb), (i7b) 

we provide numerical values of these two widths in Table I for four interesting values of 
iTi2, fn-g = 15 GeV, and sin 2 9b = 1/6. For comparison, we also present the inclusive 62 
width computed from the two-body decay matrix elements, Eq. (14). It is notable that the 
LS width is substantial in all cases, and is in fact slightly larger than the OS width for 
the lighter 62 masses we consider. The sum of the LS and OS widths, obtained from the 
three-body decay amplitudes, equals to good accuracy the inclusive width obtained from 
the two-body decay process. In general, the decay widths may depend on the sign of the 
product cos 9b sin 0&, but this dependence is absent in the limit that m& and mi vanish. 

Production of like-sign pairs, attributable directly to the Majorana nature of the gluino 
means that the subprocess of interest here generates apparent "time = zero" B° — B° flavor- 
anti-flavor mixing in e + e~ annihilation at LEP-II and linear collider energies. However, the 
&2&1 production process results in one or two jets in addition to the jets containing the b 
and b, and thus these events may not be included in a measurement of B° — B° mixing that 
focuses on bb production without additional radiation. The fraction of the decays that lead 
to like-sign pairs of 6's is 

B LS = Fls « I(9^h+3 (lg) 

r LS + r os r(s -> &i + &) + r($ &; + &) 

and from Table I, we see that this ratio is close to \ for all of the heavy bottom squark 
masses of interest, with LS slightly dominant for smaller b 2 masses. 



D. Signatures and Discovery Potential at LEP-II 

The overall process, e + e~ — > 6162, followed by 62 decay leads to a four-parton final state. 
The light bottom squarks carry color and are expected to be observed as hadronic jets. 
Absent model-dependent assumptions about bottom squark decays, these jets may have no 
special flavor content. Our SUSY subprocess results therefore in a four-jet final state: with 
2 b jets and 2 b jets. 
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The massive b\ is produced in e + e~ — > 6162 with relatively little momentum. The three 
products from its decay will therefore inherit little sense of direction. The distribution 
in cos#ij will tend to be fairly flat; subscript 1 denotes the primary b\, and j labels one 
of the decay products of the b 2 . The heavy parent Wj tosses off a b and a gluino, both 
with substantial but oppositely directed momentum. The gluino then decays into a b\ 
and the second b or b, with the daughter particles retaining the direction of the gluino's 
momentum. Since the daughter b or b follows the direction of the gluino, the two final 6's, in 
the event, whether like-sign or opposite-sign, emerge in opposite hemispheres in the overall 
e + e~ system. The invariant mass of the two final 6's will tend to be large. Furthermore, 
since rrig <C m.2, we expect the g to be highly boosted, with a small opening angle between 
its decay products. 

The two b jets are predicted to emerge in a fairly back-to-back configuration, much as 
is expected from a standard model QCD subprocess e + e~ — > (7, Z*) — > bbg, with g — > one 
or more jets. (One of the b jets may emerge fairly close to one of the b jets, resulting in 
a three-jet topology, as we investigate below.) The configuration produced by the SUSY 
process differs from that associated with e + e~ — > (7, Z*) — > qqg with g — > bb. In this later 
process, gluon splitting would yield bb pairs with modest invariant mass. 

The key question is how heavy the second scalar bottom quark 6 2 might be and still be 
discovered lurking in the LEP-II data. Alternately, one can ask what range of heavy bottom 
squark masses are ruled out by the data, if no signal is observed. We concentrate on the high 
luminosity run at y/s = 207 GeV at LEP-II, as it provides the greatest number of events for 
the masses of interest (c.f. Fig. 4). For this analysis, we do not distinguish the LS bb and 
OS bb situations, adding the distributions generated in the two cases. There would likely 
be greater potential for identifying the SUSY events if there were experimental capability 
to separate b and b jets. 

To answer our question, we address first the experimental signature of the off-diagonal 
bottom squark production process. The two bottom quark jets are almost always distin- 
guishable, with a large separation between them. Similarly, the "primary" light bottom 
squark tends to be visible as a distinct jet of hadrons. However, the light bottom squark 
from the g decay is often rather collinear with the bottom quark from the same decay, be- 
cause the g tends to be boosted by the heavy bottom squark mass. Thus, we must establish 
the number of distinctly observable jets in the final state. 
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We define an observable jet as one with transverse momentum (j>t) greater than 10 
GeV, lying in the central region of the detector, \y\ < 2, where y is the jet rapidity. The 
separation between two jets is quantified by AR = \J Ay 2 + A0 2 , where A0 is the difference 
in azimuthal angles. We consider two jets distinct from one another provided AR > 0.4; 
for smaller AR they are merged into a single jet. The distribution in the number of jets 
depends on the b 2 mass. For m 2 ~ 120 GeV, we find that the rate is split roughly evenly 
between 3 jets and 4 jets. For a heavier m 2 ~ 150 GeV, the rate is split roughly as 2/3 
3-jet events and 1/3 4-jet events. The difference arises because the larger b 2 mass results in 
a more highly boosted gluino, and thus more collinear decay products which are more likely 
to be merged into a single jet. For both masses, the 2-jet rates are smaller than a few per 
cent of the inclusive cross sections. We focus on detection of the 4-jet channel because its 
rate is a large fraction of the total rate for all masses of interest, and because we expect that 
the 4-jet configuration has smaller backgrounds, 

Backgrounds arise from e + e~ — > 3 jets and 4 jets, respectively, and involve a variety of 
mixed QCD-electroweak and purely weak processes. We simulate the 4-jet background using 
matrix elements from MADGRAPH [21]. After the acceptance cuts described above, we find 
backgrounds are typically much larger than the signal rates, hundreds of fb compared to 30 
fb (6 fb) for m 2 = 120 GeV (150 GeV). These may be reduced by positing a mass for the 
heavy bottom squark, and demanding that three of the jets reconstruct this mass within 
some window. Since the width of the heavy bottom squark is generally of order 10 GeV in 
the mass range of interest, we consider an invariant mass cut such that any three of the jets 
reconstruct an invariant mass within 10 GeV of m 2 . This m 2 -dependent cut thus forces the 
background to vary with the hypothesized value of m 2 . After its application, we find that 
the background is reduced to the manageable levels of 20 fb (43 fb) at 120 GeV (150 GeV). 

For both masses, the number of combined signal and background events is > 8, and thus 
one may apply Gaussian statistics to determine the statistical significance in the usual way, 
with a = S/y/S + B providing the confidence level (CL) for an observed signal, with S 
signal events and B background events. The resulting significances are about 5er (0.4<r) at 
m 2 = 120 GeV (150 GeV). Thus, LEP-II should be able to discover the heavier bottom 
squark through off-diagonal production if its mass is less 120 GeV. If no signal is observed, 
we estimate that masses less than 130 GeV can be excluded at the 95% CL. Our analysis 
could be improved in a number of ways, notably if experimental acceptances and efficiencies 
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were incorporated, a task beyond the scope of this work. It is our hope that the analysis 
in this paper and the exciting prospect of the discovery of SUSY will motivate a detailed 
search for signals in existing LEP-II data. 

V. CONCLUSIONS 

Mixing between weak doublets and singlets is a novel feature of the scalar fermion sector 
of the MSSM. Aside from leading to interesting phenomenology this mixing (in the top 
squark sector) is also important, allowing for large radiative corrections to the mass of the 
lightest Higgs boson, needed if the MSSM is to survive the challenge of the negative Higgs 
boson searches at LEP-II. Turning this statement around, once scalar tops are observed, 
and their masses and mixing determined, one can, to a fair degree of accuracy, determine 
whether the theory at the TeV scale is described by the MSSM, or by some more exotic 
extension. 

In this article we explore the off-diagonal squark production mode at e + e~ colliders as a 
means to measure the mixing angle and to learn more about the MSSM itself. Squark mixing 
arises from a combination of supersymmetric interactions and SUSY-breaking masses and 
trilinear terms. These in turn are related closely to the flavor structure of the MSSM. While 
there are many reasons to prefer the MSSM as the theory of physics beyond the standard 
model, its intense flavor problem indicates that SUSY-breaking is somehow very special such 
that flavor violation has not yet been observed in low energy experiments. Measurements of 
the flavor-full soft parameters will be the first experimental indications as to how nature has 
chosen to solve the SUSY flavor problem, and thus how SUSY is broken and the breaking 
communicated to the MSSM fields. 

Squark mixing is also a key player in defining the properties of the squarks themselves, 
determining the coupling to the massive electroweak bosons. Light bottom squarks, an 
interesting ingredient in the supersymmetric resolution of the bottom quark production 
cross section at hadron colliders [7], escape detection at LEP-I because their mixing angle 
is such that the left-handed and right-handed interactions with the Z boson cancel each 
other. This feature necessarily implies that off-diagonal production is non-zero and can 
be used to discover or constrain the mass of the heavier bottom squark. With a careful, 
dedicated analysis of existing LEP-II data, we show in this paper that it should be possible 
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to discover heavy bottom squarks at the 5cr level with masses as large as 120 GeV. If no 
signal is observed, exclusion limits at the 95% CL should be feasible for masses of the order 
of 130 GeV and less. Off-diagonal squark production thus allows one to explore a large 
portion of the parameter space of the light bottom squark scenario. 
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APPENDIX 

In this Appendix, we present the matrix elements for the full three-body decay subpro- 
cesses b\ — > bbb\ and b* 2 — > bbbi, keeping the dependence on the two large masses: m 2 and m~ g . 
The relevant Feynman diagram in the case of opposite sign 'OS' production (bb) is shown 
in Fig. 5 and the two diagrams for like-sign 'LS' production in Fig. 6. The 4- momenta are 
labeled as p\ and p 2 for the two bottom quarks; in the case of OS production, pi refers to 
the b and p 2 to the b. The 4-momentum of the light bottom squark is denoted p~ b . In the LS 
case, we include in \A4\ 2 the symmetry factor | for identical particles in the final state. 

The explicit expressions for these invariant amplitudes, summed/averaged over fi- 
nal/initial colors and spins are: 

\M\ls = -^j [2(l + cos 2 2^)( Pl .p g )(p 2 .p- 6 )+m?(l-cos 2 2^)( Pl .p 2 )] [| Cl | 2 + |c 2 | 2 ] 

[m~(pi • p 2 ) + 2(pi • p- b ){j>2 ■ p- b )] Re (cic* 2 ) j (19a) 



2 sin 2 29 b 



3 



8g 



M\ 2 os = ^\c 1 \ 2 {(l + cos 2 29 b )m 2 (p 1 -p 2 ) + 2(l-cos 2 2e b )(p 1 -p~ b )(p 2 -p~ b )} (19b) 
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with 



((Pi + Pi) 2 ~ m i) 2 + m l T2 ' 



(20a) 



9 9 



R r *i = ((Pi + Pi) - mpifa +Pj) ~ mp + m|T| 

6[ClC2j [((Pi + Pi) 2 - ml) 2 + mjr|] [((p 2 + pg)=» - mj)« + rojrj] ' 1 J 

To obtain the partial widths we integrate these expressions over the three body phase space 
for the 3 approximately massless final state particles. The resulting partial widths are shown 
in Table I. 
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